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摘 要 : 鼠 害 是 影响 草原 生态 健康 的 重要 因素 ,了 人 解 小 型 哮 齿 动物 种 群 密度 时 空 分 布 特征 ,对 精准 
的 鼠 害 综合 防治 具有 重要 意义 。 以 往 对 小 型 哮 齿 动物 时 空 分 布 的 研究 多 局 限于 静态 的 站 点 分 布 
或 小 范围 的 种 群 密度 时 间 变 化 分 析 , 缺 乏 对 较 大 时 空 尺 度 小 型 路 齿 动 物种 群 密度 变化 的 分 析 。 从 
已 发 表 的 文献 中 收集 了 天 山北 坡 草地 1982 一 2015 年 小 尺度 的 有 效 洞口 密度 实地 调查 信息 ,同时 结 
合 环境 因子 数据 ,再 根据 海拔 将 研究 区 划分 为 <900 m 和 >900 m2 类 ,运用 Cubist 模 型 和 随机 森林 
模型 ,分 析 有 效 洞 口 密度 时 空 分 布 。 结 果 表 明 :(1) 1982 一 2015 年 天 山北 坡 海拔 <900m 地 区 的 有 
效 洞口 密度 总 体 呈 增加 趋势 ,而 海拔 >900mm 的 地 区 总 体 呈 减少 趋势 。 基 于 Cubist 模型 构建 有 效 洞 
口 密度 与 环境 因子 的 模型 拟 合 精度 明显 优 于 随机 森林 模型 。(2) 植被 状况 、 气象 因子 和 放牧 强度 是 
天 山北 坡 有 效 洞 口 密度 时 空 分 布 主要 的 环境 驱动 因素 。 在 天 山北 坡 内 海拔 <900 m 和 >900m 的 地 
区 中 ,有 效 洞 口 密度 的 驱动 机 制 存在 着 显著 差异 。(3) 在 海拔 <900 m 地 区 ,影响 有 效 洞口 密度 时 空 


分 布 主要 是 叶 


只 指数 ,而 对 于 海拔 >900 m 地 区 为 归 一 化 植被 指数 。 这 可 能 是 受到 大 沙 鼠 (Rhom- 


bomys opimus ) 和 黄 免 尾 鼠 (Eolagurus luteus ) 消 耗 不 同类 型 植被 的 影响 。 
关 键 词 : 有 效 洞口 密度 ; 环境 因子 ; Cubist 模 型 ， 随 机 森林 模型 ， 天 山北 坡 


文章 编号 : 


哮 齿 动物 是 草地 哺乳 动物 中 种 类 和 数量 最 多 
的 一 类 动物 ,就 物质 循环 和 能 量 流动 而 言 ,其 在 草 
地 生态 系统 中 具有 重要 的 功能 和 地 位 "。 然 而 ,其 
种 群 数量 过 高 可 能 会 对 草地 造成 危害 即 草原 鼠 害 ， 
其 危害 主要 表现 在 哨 食 和 储藏 大 量 牧草 E T3 
丘 ,导致 草地 裸露 ,严重 破坏 植被 空间 结构 ,降低 草 
原生 产 力 和 草地 牧草 产量 ”一 ,打破 原 有 “ 土 - 草 - 畜 - 
BU” 的 动态 平衡 ,致使 草地 退化 严重 ”。 

新 疆 具 有 丰富 的 草地 资源 ,然而 近年 来 ,由 于 
过 度 放 牧 草地 资源 开发 不 当 全 球 气候 异常 等 ,时 
致 新 疆 草 地 鼠 害 持 续 大 面积 发 生 , 已 对 畜牧 业 的 可 
持续 发 展 和 草原 生物 多 样 性 产生 了 不 利 影响 。 为 
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了 更 及 时 地 防治 鼠 害 ,需要 加 深 对 嘴 齿 动物 种 群 密 
度 时 空 分 布 驱 动机 制 的 理解 并 更 准确 地 预测 其 时 
ZIM. HAI, BACA REI TFAA S a 
布 格局 的 研究 ,但 大 多 集中 于 人 研究 小 尺度 的 种 群 数 
量 随时 间 的 变化 ”或 研究 其 潜在 空间 分 布 , 却 很 少 
涉及 大 尺度 的 种 群 密度 时 空 分 布 的 变化 特征 及 其 
机 制 。 如 乔 雪 丽 等 "利用 最 大 焙 模 型 (MaxFEnb) 对 在 
不 同 环境 条 件 下 的 大 沙 鼠 在 亚洲 的 潜在 地 理 分 布 
进行 预测 ; 何 咏 琪 等 "基于 "3S ”技术 与 地 面 实地 调 
查 数 据 , 建 立 了 草原 鼠 害 监 测 模型 并 模拟 出 了 青海 
省 草原 鼠 害 分 布 区 ;Mohammadi 等 中 模拟 了 伊朗 气 
修 变 化 下 两 种 沙漠 跳 鼠 的 当前 和 未 来 潜在 分 布 。 
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E 态 效应 .干旱 区 主要 生态 系统 对 气候 变化 和 人 类 活动 的 响应 模拟 等 
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TR WES. 基于 Cubist 模 型 的 天 山北 坡 草地 鼠 群 密度 时 空 分 布 特征 


目前 ,环境 因子 已 被 广泛 用 于 物种 种 群 密度 分 
布 的 驱动 机 制 分 析 和 时 空 分 布 预测 ”“"。 对 距 齿 动 
物 而 言 ,环境 条 件 在 其 种 群 密度 波动 和 空间 分 布 中 
扮演 着 重要 角色 ,可 通过 直接 或 其 相互 作用 对 其 产 
生 重 要 影响 ”…。 嘴 雌 动 物种 群 密度 的 时 空 分 布 受 地 
形 中 土壤 质 地 下 气温 ”降水 “植被 "放牧 
等 的 影响 ,上 且 其 具有 年 际 和 季节 性 变化 的 特征 ”。 
因此 ,通过 建立 路 此 动物 种 群 密度 与 环境 因子 的 关 
系 来 预测 其 种 群 密度 的 时 空 分 布 具 有 一 定 可 行 
性 。 在 衡量 害 鼠 种 群 密度 的 指标 中 ,有 效 洞 口 密度 
是 最 有 效 的 指标 之 一 "”"。 该 指标 已 广泛 地 在 实地 
调查 中 被 估计 并 记录 ,样本 数 相对 丰富 ,适用 于 本 
研究 中 大 尺度 种 群 密度 的 归 一 性 分 析 和 预测 。 

天 山 山脉 北 坡 属于 中 亚 干 旱地 区 上 典型 的 山地 - 
绿洲 -沙漠 (从 高 地 到 低地 ) 生 态 系统 , 动 植物 的 
分 布 很 大 程度 上 取决 于 海拔 的 变化 。 在 天 山北 坡 ， 
嘴 齿 动物 对 草原 的 损害 主要 由 分 布 在 平原 荒漠 草 
地 的 大 沙 鼠 (Rhombomys opimus ) 和 山地 草原 的 黄 免 
FÆ M (Eolagurus luteus) 引 起 , 旦 大 沙 鼠 主要 栖息 在 
EIR 900 m LA P EU HEHE, ALI, AS CARA XC 
献 中 数据 特点 将 天 山北 坡 分 为 海拔 <900 m 和 >900 
m2 类 区 域 , 且 主 要 的 优势 种 分 别 为 大 沙 鼠 和 黄 锡 
FÉ BR S 

25 E. ,本文 以 有 效 洞口 密度 反映 害 鼠 的 种 群 密 
度 , 从 已 发 表 的 文献 收集 天 山北 坡 草地 小 尺度 的 有 


阿拉 湖 he. "S 


Af ils] H RE BE RI S H8] As ROR SEIT 
分 析 等 手段 ,对 2 类 区 域 分 别 运 用 随机 森林 (RF) 模 
型 和 Cubist 模 型 模拟 有 效 洞口 密度 。 最 后 ,选取 较 
优 方法 定量 研究 环境 因子 对 有 效 洞口 密度 时 空 分 
布 的 影响 ,对 1982 一 2015 年 的 天 山北 坡 草地 有 效 洞 
口 密度 时 空 分 布 进行 预测 ,并 分 析 其 变化 特征 。 


1 数据 与 方法 


1.1 研究 区 概况 

天 山北 坡 ( 图 1) 是 指 由 横贯 新 疆 中 部 天 山 山 肪 
发 育 的 河流 向 北 流 经 的 区 域 , 地 理 位 置 为 79"53'~ 
96°06'EE,42°50'~46°12'N, 东 起 伊 吾 , 西 至 温泉 县 ， 
北 连 古 尔 班 通 古 特 沙漠 。 区 域内 自然 条 件 的 垂直 
分 带 性 明显 ,地势 南 高 北 低 ,海拔 由 山地 的 5000 m 
左右 下 降 到 沙漠 边缘 的 200 m 左右。 由 于 位 于 亚 欧 
大 陆 腹地 ,远离 海洋 ,属于 典型 的 温带 大 陆 性 气候 ， 
其 特点 是 冬季 寒冷 漫长 ,夏季 干旱 炎热 ,年 均 温 在 
-6.5-11.0 % ,降水 季节 变化 显著 ,年 内 降水 量 主要 
集中 在 5.6 月 ,2 月 最 少 ,年 降水 量 由 南部 山区 中 低 
山 带 的 500 mm 降 至 北部 沙漠 区 的 100 mm 左右 。 天 
山北 坡 独 特 的 自然 条 件 使 其 成 为 全 球 温带 干旱 区 
大 型 山地 -绿洲 -沙漠 生态 系统 的 典型 代表 ?2 。 
1.2 数据 来 源 

本 研究 采用 的 数据 主要 包括 有 效 洞 口 密度 和 


图 1 基于 文献 搜集 的 天 山北 坡 鼠 洞 调 查 样 点 分 布 


Fig. 1 Sample distribution of rat hole investigation on the north slope of Tianshan Mountains based on literature collection 
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环境 因子 数据 ( 表 1)。 有 效 洞口 密度 数据 通过 已 公 
开发 表 的 文献 汪汪 获取 ,所 搜集 记录 中 仪 有 地 名 但 
无 经 纬度 信息 的 通过 GPSSPG 网 站 (http://www. 
gpsspg.com/) 进 行 查询 ”“'。 对 得 到 的 数据 进行 整理 ， 
最 终 得 到 了 275 条 具有 时 空 信息 的 有 效 洞口 密度 数 
据 ,年 份 跨度 为 1982 一 2015 年 。 本 文 初步 选择 了 21 


1.3 研究 方法 

首先 ,通过 ArcGIS 10.2 软件 对 获得 的 1982— 
2015 年 的 有 效 洞口 密度 时 空 信息 进行 整理 ,生成 含 
有 准确 空间 位 置信 息 的 月 尺度 矢量 文件 。 随 后 , 根 
据 每 个 样本 的 时 空 信息 提取 对 应 的 21 个 环境 因子 
并 根据 海拔 将 天 山北 坡 分 为 <900 m 和 >900 m 2 


个 可 能 影响 有 效 洞口 密度 的 环境 因子 ,主要 包括 年 类。 最 后 ,分 别 基 于 RF 模型 和 Cubist 模 型 拟 合 有 效 


43 月份. 地 形 .土壤 质地 植被 .气温 、 降 水量 ,短波 


洞口 密度 ,并 选取 较 优 模型 定量 研究 有 效 洞口 密度 


辐射 和 放牧 强度 。 时 空 分 布 的 变化 特征 及 其 机 制 。 
表 1 环境 因子 数据 来 源 
Tab. 1 Data sources of environmental factors 
环境 因子 空间 分 辩 率 数据 来 源 

海拔 30 m 美国 地 质 勘探 局 全 球 卫星 影像 (USGS Earth Explorer, http://earthexplorer.usgs.gov ) 

坡度 30 m 美国 地 质 勘 探 局 全 球 卫 星 影像 (USGS Earth Explorer, http://earthexplorer.usgs.gov ) 

坡 向 30 m 美国 地 质 勘 探 局 全 球 卫 星 影像 (USGS Earth Explorer, http://earthexplorer.usgs.gov ) 

粉 砂 含量 ] km 国家 冰川 冻 土 沙漠 科学 数据 中 心 (http:/www.ncdc.ac.enjm 

WERE 1 km 国家 冰川 冻 土 沙漠 科学 数据 中 心 (http:/www.ncde.ac.emm 

黏土 含量 ] km 国家 冰川 冻 土 沙漠 科学 数据 中 心 (http:/www.ncdc.ac.enjm 

归 一 化 植被 指数 (NDVI1) 8 km 美国 国家 航天 航空 局 (NASA ) 的 全 球 监测 与 模型 研究 组 (GIMMS ) 提 供 的 3g.v1 数 
据 (https://ecocast.arc.nasa.gov/data/pub/gimms/) 

叶 面 积 指数 (LAI) 8 km NOAA 人 气候 数据 记录 (CDR ) 的 AVHRR 叶 面 积 指数 (LAIT) 和 吸收 光合 有 效 辐 射 的 
部 分 (FAPAR ) 数 据 集 

年 平均 气温 1 km 中 国 1 km 分 辩 率 逐 月 平均 气温 数据 集 (1901 一 2017 年 )” 

年 最 低 气 温 1 km 中 国 1 km 分 辨 率 月 最 低温 度数 据 集 (1901 一 2017 年 )*™ 

年 最 高 气温 1 km 中 国 1 km 分 辩 率 月 最 高 温度 数据 集 (1901 一 2017 年 )*" 

月 平均 气温 1 km 中 国 1 km 分辨 率 逐 月 平均 气温 数据 集 (1901 一 2017 年 ) 

月 最 高 气温 1 km 中 国 1 km 分 辩 率 月 最 高 温度 数据 集 (1901 一 2017 年 )*" 

月 最 低 气温 1 km 中 国 1 km 分 辩 挛 月 最 低温 度数 据 集 (1901 一 2017 年 ) 

年 累积 降水 量 1 km 中 国 1 km SP SERIE KEk BOSE SE (1901—2017 4F 

调查 前 3 月 的 平均 降水 量 1 km 中 国 1 km 分 辨 率 逐 月 降水 量 数据 集 (1901 一 2017 年 )%” 

月 累积 降水 量 1 km 中 国 1 km 42358248 HE KEk HESS (1901—2017 4) 

放牧 强度 0.08333° 世界 粮农 组 织 (FAO, http://www.fao.org/livestock-systems/en/) 野 外 实测 数据 ,新疆 
appe 

短波 辐射 0.1° 中 国 区 域 地 面 气象 要 素 驱动 数据 集 (1979 一 2018 年 ) 

1.3.1 REF 模型 REF 模型 是 一 种 以 决策 树 为 基本 单 。” 格 搜索 算法 "探索 各 区 数据 的 最 优 参 数组 合 ,采用 


元 的 集成 学 习 方法 “。 其 基于 Bagging 算 法 ,从 原始 
样本 集 N 中 有 放 回 地 随机 抽取 n 个 样本 生成 新 的 训 
练 样本 集 ,通过 算法 训练 得 到 "” 棵 决策 树 ,最 后 通过 
平均 的 方式 整合 各 决策 树 的 预测 结果 , 便 可 得 到 RF 
模型 的 最 终结 果 。RF 模 型 在 每 个 回归 树 的 建立 过 
程 中 ,对 样本 空间 和 特征 空间 进行 随机 采样 ,引入 
随机 属性 ,进而 降低 了 回归 树 模 型 之 间 的 相关 性 ， 
又 通过 组 合 大 量 的 回归 树 ,改善 了 模型 的 泛 化 性 ， 
使 得 此 算法 具备 高 效 、 准 确 等 特点 。 本 研究 利用 网 


Python 中 scikit-learn 包 进行 模型 构建 '*。 

1.3.2 Cubist 模 型 ”Cubist 模 型 是 Quinlan 的 M5' 模 
型 树 的 进一步 发 展 »。Cubist 模 型 在 多 元 线性 模型 
中 得 到 一 系列 让 then 形 式 的 规则 , 当 一 组 变量 符合 
规则 的 条 件 时 ,使 用 相应 的 模型 计算 变量 的 预测 
值 "'。Cupbist 算 法 用 一 系列 组 合 起 来 的 分 段 线性 模 
型 ,很 好 地 解决 了 非 线 性 问题 。 此 外 ,该 方法 还 具 
有 变量 的 重要 性 分 析 功 能 “'。 本 研究 使 用 R 软 件 中 
Cubist 程序 包 ,其 主要 参数 为 committees 和 neigh- 
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13.3 模型 精度 验证 ”本文 从 数据 集 随 机 选择 80% 
作为 模型 的 训练 样本 集 ,剩余 的 20% 作 为 模型 验证 
样本 和 集 。 选 用 均 方 根 误差 (RMSE) 和 决定 系数 (R) 
作为 模型 精度 评价 指标 ,综合 评价 RF 模型 和 Cubist 
模型 的 预测 能 力 。 计 算 公式 如 下 : 


RMSE = (1) 


gr-Ym-Yy/Ymy-1y 2 


SUSY, Y; 和 了 分 别 为 有 效 洞口 密度 实际 值 预 
测 值 和 样本 平均 值 ; n 为 样本 数 。 


2 结果 与 分 析 


2.1 有 效 洞口 密度 样本 数据 分 布 特征 分 析 

对 从 文献 中 获取 的 有 效 洞 口 密度 数据 进行 统 
计 分 析 , 如 表 2 所 示 , 共 275 条 数据 ,其 中 最 小 值 为 1 
个 .hm ,最 大 值 为 1370 个 .hm ,平均 值 为 254.23 
个 .hm?。 海拔 <900 m 地 区 的 样 点 有 效 洞口 密度 范 
围 为 4~1220 个 .hm ,海拔 >900 m 地 区 的 样 点 有 效 
洞口 密度 在 1~1370 个 .hm-? 之 间 。 
2.2 草原 有 效 洞口 密度 拟 合 模型 精度 评价 与 优选 

将 全 部 环境 因子 作为 自 变量 ,有 效 洞口 密度 为 
VS AE t ,分别 构建 RF 模 型 和 Cubist 模 型 并 计算 尼 和 
RMSE( 图 2)。 结 果 表 明 ,在 海拔 <900m 地 区 其 RF 
模型 和 Cubist fo 78 HY R^ 3 3] Jy 0.773 和 0.832 ,在 海 
£ 2900 m Hi X F; RF $R AAN Cubist 模型 的 尼 均 在 
0.600 以 上 , Cubist 模型 均 高 于 RF 模型 。 同 时 ,RF 模 
型 的 RMSE 均 大 于 Cubist 模 型 。 综 上 ,Cubist 模 型 具 
有 更 高 的 拟 合 精度 和 更 小 的 误差 ,因此 本 研究 采用 
Cubist 模型 来 预测 天 山北 坡 草地 有 效 洞口 密度 。 
2.3 影响 因子 的 重要 性 分 析 

本 文 基于 Cubist 模 型 各 因子 的 重要 性 排序 ,分 
析 环 境 因 子 对 天 山北 坡 草地 有 效 洞口 密度 时 空 分 
布 的 影响 (图 3)。Cubist 模 型 可 根据 样本 信息 将 相 
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ATE AS SR E198 A KFAR ,结合 预测 因子 的 相对 重 
要 性 ,筛选 出 预测 结果 较 好 的 影响 因子 。 重 要 性 显 
示 海 拔 < 900 m 和 海拔 >900 m 地 区 之 间 存 在 着 显著 
差异 。 在 天 山北 坡 海拔 <900 m tl DC ,对 草地 有 效 洞 
口 密度 分 布 贡献 最 大 的 是 月 最 低 气 温 , 其 次 为 月 平 
均 气 温 、 叶 面积 指数 (LAIT)。 而 对 于 >900 m 地 区 , 贡 
献 最 大 的 为 年 累积 降水 量 , 其 次 为 放牧 强度 、 归 一 
化 植被 指数 (NDVI)。 气 候 因 子 .植被 状况 与 人 类 活 
动因 子 对 天 山北 坡 有 效 洞口 密度 时 空 分 布 影响 较 大 。 
2.4 天 山北 坡 草 地 有 效 洞口 密度 时 空 分 布 

使 用 此 Cubist 模 型 对 研究 区 草地 有 效 洞口 密度 
时 空 分 布 进行 反 演 (图 4)。 并 根据 Cubist 模 型 分 别 
计算 1982 一 2015 年 天 山北 坡 草 地 有 效 洞口 密度 在 4 
月 和 10 月 的 平均 值 (图 5)。 比 较 1982 .1993 .2004 
年 和 2015 年 有 效 洞 口 密度 的 年 际 空间 分 布 差异 以 
及 其 春季 (4 月 ) 和 秋季 (10 月 ) 的 季节 性 变化 。 结 
表明 ,有 效 洞口 密度 具有 明显 的 区 域 差 异 ,靠近 人 
类 活动 频繁 的 农田 生态 系统 周围 草地 有 效 洞口 密 
度 普 遍 较 高 ,而 北部 受 人 类 干扰 较 少 的 车 江 地 区 其 
密度 较 低 。 在 海拔 <900m 地 区 4 月 有 效 洞口 密度 
平均 为 354 个 .hm ,10 月 平均 为 424 个 :hm2 ;海拔 > 
900 m 地 区 4 月 有 效 洞 口 密度 平均 为 421 个 hm”,10 
月 平均 为 407 个 hm”*。1982 一 2015 年 海拔 < 900 m 
地 区 有 效 洞口 密度 整体 为 波动 递增 的 趋势 ,34a 间 4 
月 的 有 效 洞 口 密度 增加 了 37.25% , 10 月 增加 了 
41.5896 ; 而 海拔 >900 m 地 区 整体 为 波动 递减 的 趋 
势 ,34a 间 4 月 减少 了 4.83% ,10 月 减少 了 26.99% , 
且 草 地 有 效 洞口 密度 波动 剧烈 (图 5)。 在 天 山北 坡 
海拔 <900 mm 地 区 每 年 10 月 的 有 效 洞口 密度 相对 于 
4 月 显著 增加 ,体现 出 较 强 的 季节 差异 性 。 


3 讨论 


本 研究 主要 基于 Cubist 模 型 .结合 环境 因子 预 
测 有 效 洞口 密度 ,基本 实现 了 草原 小 型 路 齿 动 物种 
群 密度 时 空 量化 与 表达 ,有 效 地 弥补 了 目前 在 大 时 


表 2 有 效 洞口 密度 统计 描述 


Tab. 2 Statistical description of active burrow entrance densities 


数据 类 型 样本 数 最 小 值 /个 hm? 
全 部 数据 275 1 
海拔 <900 m 地 区 数据 190 4 
海拔 >900 m 地 区 数据 85 1 


最 大 值 /个 .hm” 平均 值 /个 hm” 标准 差 
1370 254.23 254.23 
1220 251.76 222.54 
1370 259.73 253.62 
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已 有 的 研究 因 缺 乏 长 时 间 序 列 以 及 较 大 研究 区 的 
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图 2 有 效 洞口 密度 拟 合 值 在 海拔 <900 m 和 >900 m 地 区 的 验证 散 点 图 
Fig.2 Scatter polt of fitting values of active burrow entrance densities in areas with an elevation 
less than or equal to 900 m and higher than 900 m 
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注 :LAI 为 叶 面积 指数 ;NDVI 为 归 一 化 植被 指数 。 
图 3 有 效 洞口 密度 的 影响 因子 重要 性 
Fig. 3 Importance of influencing factors for active burrow entrance densities 
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的 认识 存在 局 限 性 
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4 1982 一 2015 年 天 山北 坡 草 地 每 11 a 的 有 效 洞口 密度 空间 分 布 
Fig.4 Spatial distributions of active burrow entrance densities in grassland on the north slope of Tianshan Mountains 
every 11 years from 1982 to 2015 
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图 5 1982 一 2015$ 年 天 山北 坡 草地 有 效 洞 口 密度 在 海拔 < 900 m 和 >900 m 地 区 的 变化 趋势 


Fig. 5 Change trend of active burrow entrance densities in grassland on the north slope of Tianshan Mountains in areas with an 


elevation less than or equal to 900 m and higher than 900 m from 1982 to 2015 
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在 研究 方法 方面 ,已 有 的 研究 中 鲜 有 将 RF、 
Cubist 模 型 等 机 器 学 习 方 法 应 用 到 小 型 路 肯 动 物种 
群 密度 估算 的 研究 方面 。 因 此 ,本 研究 尝试 将 不 同 
机 需 学 习 方 法 在 该 领域 进行 应 用 ,这 也 是 本 文 的 主 
要 创新 点 之 一 。 所 构建 的 Cubist 模 型 在 量化 环境 因 
子 对 有 效 洞口 密度 的 影响 方面 , 优 于 REF 模 型 。 草 
原 小 型 嘴 此 动 物 生 态 特 点 是 种 群 数量 波动 剧烈 ,有 
些 年 份 较 低 ,有些 年 份 猛 增 , 且 其 分 布 具 有 明显 的 
区 域 性 差异 2 。Cubist 模型 可 以 针对 不 同 区 域 构建 
合适 的 应 用 模型 ,其 能 够 同时 处 理 连 续 的 和 离散 的 
数据 ,因此 ,相对 于 RF 模型 ,Cubist 模 型 可 能 能 
更 好 地 模拟 草原 有 效 洞 口 密度 。 

首先 ,本 研究 发 现 植被 状况 在 鼠 类 种 群 密度 的 
时 空 分 布 中 占 重要 地 位 ,但 一 个 有 趣 的 现象 是 ,在 
海拔 < 900 m 地 区 ,影响 有 效 洞口 密度 时 空 分 布 的 主 
要 是 LAI, 而 对 于 海拔 >900 m 地 区 为 NDVI( 图 3)。 
这 可 能 是 由 于 遥感 NDVI 数 据 在 监测 梭 梭 的 生长 和 
地 下 生物 量 方面 不 是 很 有 效 。 因 为 在 海拔 <900 m 
的 平原 荒漠 地 区 ,大 沙 鼠 将 梭 梭 的 嫩 根 和 树枝 作为 
主要 食物 ,而 黄 兔 尾 鼠 主要 食用 青 某 属 植物 。 有 效 
洞口 密度 与 LAI 或 NDVI 之 间 关 系 的 差异 可 能 受到 
海拔 <900 m 和 >900 mm 地 区 中 大 沙 鼠 和 黄 兔 尾 鼠 所 
消耗 不 同类 型 植被 的 影响 。 植 被 生长 状况 直接 决 
定 鼠 类 食物 数量 和 品质 中。 植物 群落 的 各 项 指标 
(植被 高 度 、 盖 度 、 密 度 和 地 下 生物 量 以 及 群落 多 样 
性 等 ) 都 综合 地 影响 着 中 上 雳 动物 种 群 密度 “” S DS 
本 文中 有 效 洞口 密度 数据 是 通过 已 公开 发 表 的 文 
献 获 取 ,而 文献 中 缺乏 调查 有 效 洞 口 密度 时 所 涉及 
的 植被 基本 性 状 指 标 , 且 我 们 目前 无 法 获取 天 山北 
坡 植被 群落 的 空间 数据 。 因 此 ,本文 环境 变量 中 的 
植被 性 状 指 标 只 涉及 植被 指数 (NDVI、LAI) , 虽 植 被 
指数 已 广泛 用 来 定性 和 定量 评价 植被 覆盖 及 其 生 
长 活力 ,但 对 于 优势 植物 种 类 等 性 状 指标 的 缺 
= ,可 能 会 使 环境 因子 的 预测 能 力 受 到 很 大 的 限 
制 。 此 外 ,使 用 NDVI 和 LAI 数 据 较 粗糙 的 空间 分 辩 
率 可 能 会 给 此 类 分 析 华 来 不 确定 性 ,特别 是 在 植被 
和 覆盖 率 较 低 的 沙漠 地 区 ,这 可 能 导致 低估 空间 植被 
变异 性 对 有 效 洞 口 密度 的 影响 。 除 食物 因素 外 , 气 
候 因 子 与 放牧 强度 对 嘴 齿 动物 种 群 密度 也 产生 了 
极 大 的 影响 ,并 且 影 响 会 因 地 区 而 异 。 对 比 2 类 地 
区 ,发 现 有 效 洞口 密度 对 温度 响应 显著 的 区 域 为 海 
J< 900 m 地 区 ,而 降水 影响 显著 的 区 域 为 >900 m 


地 区 。 这 可 能 是 因为 在 >900 m 海 拔 带 降水 的 时 空 
异 质 性 较 高 ,因此 对 植被 生长 的 影响 更 大 ', 这 与 
天 山北 坡 和 全 国 范 围 的 相关 性 空间 人 研究 结论 是 相 
吻合 的 "。 此 外 , 梭 梭 是 大 沙 鼠 的 主要 食物 来 源 ， 
而 黄 免 尾 鼠 主要 食用 青 蔽 属 植物 , 且 黄 兔 尾 鼠 其 在 
阴雨 刊 风 等 气温 较 低 天 气 条 件 下 ,活动 性 明显 降 
低 。 这 可 能 也 是 黄 兔 尾 鼠 种 群 密度 对 降水 更 敏感 
的 原因 。 在 自然 情况 下 , 鼠 类 通过 与 草地 环境 的 相 
互 作用 、 相 互 协调 使 之 经 常 保持 动态 的 平衡 状态 。 
然而 ,长 期 以 来 ,人 们 为 经 济 利益 所 驱动 的 高 强度 
利用 方式 ,深刻 地 影响 鼠 类 及 天 敌 动 物 赖 以 生存 地 
草地 基质 条 件 ,而 群居 的 植 食性 鼠 类 生殖 力 强 ,种 
群 数量 增长 快 ,往往 在 此 类 次 生 环 境 中 占据 优势 而 
发 生 鼠 害 允 。 一 些 学 者 基于 大 面积 的 草地 调查 资 
料 也 发 现 了 草地 鼠 类 种 群 密度 的 空间 分 布 与 放牧 
TREE A RA 

其 次 ,天 山北 坡 海拔 >900 m 靠 近 人 类 活动 频繁 
的 农田 生态 系统 周围 的 草地 区 域 ,其 有 效 洞口 密度 
普遍 较 高 ,这 可 能 是 因为 黄 免 尾 鼠 主要 栖息 在 丘陵 
及 总 漠 草 原 ,农垦 后 聚集 在 农田 及 水 渠 附近 '”。 黄 
兔 尾 鼠 数量 年 际 变动 较为 剧烈 ,往往 在 中 低 密 度 水 
平 连续 维持 若干 年 之 后 ,种群 密度 突然 升 高 ,而 后 
又 处 于 低潮 期 。 低 潮 后 ,在 自然 条 件 适宜 的 年 份 密 
度 又 开始 增 大 (图 5)。 由 本 研究 可 以 发 现 , 黄 兔 尾 
鼠 数量 的 年 际 变动 可 能 主要 归 因 于 降水 量 .放牧 强 
度 和 植被 状况 ,但 鉴于 黄 兔 尾 鼠 种 群 密度 与 植被 - 
降水 -人 类 活动 的 关系 复杂 、 时 空 不 确定 性 大 ,这 一 
现象 的 产生 机 制 ,还 需要 进一步 探讨 。 

最 后 ,环境 因子 和 有 效 洞口 密度 是 一 个 动态 变 
化 的 过 程 ,遥感 技术 能 够 及 时 进行 大 范围 监测 了 解 
环境 各 个 时 期 的 变化 。 有 效 洞口 密度 的 时 空 动 态 
在 响应 环境 因子 和 人 类 活动 的 复杂 作用 时 展现 出 
ADMIRU ,需要 从 多 视角 研究 小 型 咕 齿 动物 种 群 密 
度 时 空 变化 。 本 文选 择 了 1982 一 2015 年 的 数据 进 
行 建 模 , 能 够 反映 路 齿 动 物种 群 密度 的 年 际 和 季节 
的 变化 情况 ,但 对 于 偶尔 极 高 的 繁殖 率 ,环境 因子 
的 预测 能 力 可 能 受到 很 大 的 限制 。 此 外 , 鼠 类 种 群 
年 龄 结构 .天 敌 种 群 分 布 以 及 人 为 因素 “7 在 鼠 洞 
密度 分 布 的 预测 中 具有 重要 影响 ,这 在 未 来 的 研究 
中 需要 更 多 关注 。 


4 结论 


本 研究 收集 了 1982 一 2015 年 的 天 山北 坡 有 效 
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洞口 密度 的 小 尺度 时 空 分 布 数据 和 环境 因子 数据 ， 
并 基于 Cubist 模型 估算 天 山北 坡 草地 有 效 洞口 密 
BE ,量化 了 环境 要 素 对 小 型 中 此 动物 种 群 密度 的 影 
啊 ,取得 了 较 好 的 研究 精度 和 效果 ,为 区 域 草原 鼠 
害 预防 预测 工作 提供 了 新 思路 .新 方法 。 得 到 以 下 
结论 : 

(1) 利用 Cubist 模 型 和 RF 模型 对 天 山北 坡 草 
地 有 效 洞 口 密度 进行 预测 ,海拔 < 900 m 地 区 拟 合 效 
果 尼 分 别 为 0.832 和 0.773 ,海拔 >900 m HX RAYI 
为 0.667 和 0.612, Cubist 模型 拟 合 精度 明显 优 于 RF 
模型 。 

(2) 基于 Cubist 模 型 的 影响 因子 重要 性 分 析 结 
果 ,植被 状况 .气象 因子 和 放牧 强度 是 天 山北 坡 有 
效 洞 口 密度 时 空 分 布 的 主要 环境 驱动 因素 。 在 天 
山北 坡 内 的 海拔 <900 m 和 >900 m 地 区 中 ,有 效 洞 
口 密度 的 驱动 机 制 存在 着 显著 差异 。 

(3) 在 海拔 < 900 m 地 区 ,气温 对 有 效 洞口 密度 
的 贡献 度 比 在 海拔 >900 m 地 区 更 高 ,而 降水 量 对 其 
影响 更 低 。 在 海拔 < 900 m 地 区 ,影响 有 效 洞口 密度 

时 空 分 布 的 主要 是 LAI, 而 对 于 海拔 >900 m 地 区 为 
NDVI。 这 可 能 是 受到 大 沙 鼠 和 黄 兔 尾 鼠 所 消耗 不 
同类 型 植被 的 影响 。 

(4) 有 效 洞口 密度 的 时 空 分 布 具有 明显 的 年 份 
和 季节 波动 ,1982 一 2015 年 海拔 <900 m 地 区 有 效 
洞口 密度 整体 为 波动 递增 的 趋势 ,34a 间 4 月 有 效 
洞口 密度 增加 了 37.25% , 10 月 增加 了 41.58% ; 而 海 
4/2900 mm 地 区 整体 为 波动 递减 的 趋势 ,34a 间 4 月 
减少 了 4.83% ,10 月 减少 了 26.99%。 

(5) 在 海拔 >900 m 地 区 ,草地 有 效 洞口 密度 波 
动 剧 烈 , 有 效 洞口 高 密度 区 主要 分 布 在 靠近 人 类 活 
动 的 牧区 和 绿洲 区 。 这 可 能 主要 归 因 于 降水 量 、 放 
牧 强度 和 植被 状况 的 著 加 效应 ,但 鉴于 黄 兔 尾 鼠 种 
群 密度 与 植被 -降水 -人 类 活动 的 关系 复杂 时空 不 
确定 性 大 ,这 一 现象 的 产生 机 制 , 还 需要 进一步 探 
讨 。 未 来 可 通过 大 数据 分 析 手 段 提高 数据 的 精度 
或 多 视角 考虑 种 群 密度 的 影响 因子 ,尤其 是 植被 基 
本 性 状 指标 ,并 与 空间 异 质 性 探测 方法 结合 ,多 层 
次 深化 本 研究 方向 和 主题 
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Spatiotemporal distribution of rat population density in grassland on the 
north slope of Tianshan Mountains based on Cubist model 


PAN Qun", SHI Haiyang”, ZHANG Wenqiang", LUO Geping"" CHEN Chunbo^ 
(1. College of Resources and Environmental Sciences, Xinjiang University, Urumqi 830046, Xinjiang, China; 2. State Key 

Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi 
830011, Xinjiang, China; 3. University of Chinese Academy of Sciences, Beijing 100049, China; 4. Central Asia Center for 


Ecology and Environmental Research Center, Chinese Academy of Sciences, Urumqi 830011, Xinjiang, China) 


Abstract: Rodent damage is an important factor affecting grassland ecological health. Understanding the 
spatiotemporal distribution characteristics of small rodent population density is crucial for accurate 
comprehensive rodent control. Previous studies on the spatiotemporal distribution of small rodents were mostly 
limited to static site distribution or time change analysis of small-scale population density; investigation on the 
population density change of small rodents on a large temporal and spatial scale is lacking. In this study, the field 
survey information of small-scale active burrow entrance densities on the grassland on the north slope of Tianshan 
Mountains from 1982 to 2015 was collected from the literature and combined with environmental variable data. 
The study area was divided as areas with elevation of <900 m and >900 m. The spatiotemporal distribution of 
active burrow entrance densities was analyzed using the Cubist and random forest (RF) models. The following 
results were obtained. (1) From 1982 to 2015, the active burrow entrance densities increased in the area with an 
elevation of <900 m on the north slope of Tianshan Mountains, whereas that in the area with an elevation of >900 
m decreased. The fitting accuracy of the Cubist model of active burrow entrance densities and environmental 
factors was better than that of the RF model. (2) Vegetation status, meteorological factors, and grazing intensity 
were found to be the primary environmental driving factors for the spatiotemporal distribution of active burrow 
entrance densities on the north slope of Tianshan Mountains. Significant differences in the driving mechanism of 
active burrow entrance densities were found between the areas with an elevation of <900 m and those with an 
elevation of >900 m on the north slope of Tianshan Mountains. (3) The main factor affecting the spatiotemporal 
distribution of active burrow entrance densities was the leaf area index and normalized difference vegetation 
index for the areas with elevation <900 m and >900 m, respectively. This difference may be affected by the 
consumption of different types of vegetation by gerbils (Rhombomys opimus) and yellow rabbit tail rats 
(Eolagurus luteus). 

Key words: active burrow entrance densities; environmental factors; Cubist model; random forest model; 


northern slope of Tianshan Mountains 


